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Abstract 

METATOYs can change the direction of light in ways that appear to, but do not 
actually, contravene the laws of wave optics. This direction change applies only to 
part of the transmitted light beam; the remainder gets re-directed differently. For a 
specific example, namely confocal pairs of rectangular lenslet arrays with no dead 
area between lenslets, we calculate here the fractions of power of a uniform-intensity 
light beam incident from a specific (but arbitrary) direction that get re-directed 
in different ways, and we derive an equation describing this redirection. This will 
facilitate assessment of the suitability of METATOYs for applications such as solar 
concentration. Finally, we discuss similarities between the multiple refraction of 
light at the lenslet arrays and multiple refraction and reflection of cold atoms at a 
barrier in the presence of the light fields. 

Key words: confocal lenslet arrays, METATOYs, field of view, geometrical optics, 
optical materials 
PACS: 42.15.-i 



1 Introduction 

METATOYs pQ can be described as windows that "refract" (change the direc- 
tion of) transmitted light rays. This direction change can be, for example, a 
rotation around the window normal [2J, or indeed around any other direction 
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Fig. 1. Confocal lenslet arrays (CLAs). The focal lengths of the lenslets in the left 
and right array are respectively fi and fi- (a) Example of a light ray passing through 
corresponding lenslets (solid arrow) and non-corresponding lenslets (dashed arrow). 
a and f3 are the angles of incidence and refraction, respectively, of the former light 
ray. (b) Imaging properties of light rays that pass through corresponding lenslets in 
planar CLAs. 

[3] ; a flipping (sign change) of one of the transverse direction components [I] ; 
negative refraction that leads to pseudoscopic imaging [SJ; and a variation of 
Snell's law in which sines are replaced by tangents [B] (and a generalisation 
thereof [7]). It can be shown that most of these light-ray-direction changes 
can lead to light-ray fields that cannot be represented wave-optically [Tf8] . 
However, this is not actually the case, and this apparent conflict with wave 
optics is resolved by introducing discontinuities into the wave front [9"j. 

The wave-front discontinuities introduced by METATOYs result in only part 
of the light being redirected as advertised; the remainder undergoes a different 
direction change. This imperfection can be remedied by more careful optical 
design, for example insertion of arrays of field lenses into the common focal 
plane of confocal lenslet arrays (CLAs) pU], which are examples of META- 
TOYs j6]. Nevertheless, it compromises the performance of METATOYs as 
components in visual optical instruments and in instruments for light-shaping 
applications. 

One timely application of light shaping is solar concentration [11]. Testing 
whether or not METATOYs can offer anything new in this important area 
is an obvious task, one for which it is crucial to understand how much light 
undergoes what direction change. 

Here we address this point for simple CLAs [6] (Fig. [I]). Note that lenslet (or 
microlens) arrays are already being investigated as a way of improving the 
efficiency of photovoltaic cells (e.g. [12] )• We build on previous work on the 
mechanism behind, and qualitative effect of, different parts of a light beam 
being refracted differently. Specifically, we calculate the fraction of a uniform- 
intensity light beam, incident from a specific (but arbitrary) direction, that un- 
dergoes the "correct" direction change upon transmission through the CLAs; 
we call this fraction (. We also calculate the fractions ( m of transmitted light 
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Fig. 2. Geometries of CLAs to which our calculation applies, (a) Confocal arrays of 
cylindrical lenslets; (b) confocal rectangular arrays of spherical lenslets with rectan- 
gular apertures. In both cases, f\ and fi is the focal length of the lenslets in the back 
array and the front array, respectively, and there is no dead area between neigh- 
bouring lenslets. For clarity, a pair of lenslets that occupy corresponding positions 
in the two lenslet arrays are highlighted in both diagrams. 



that undergo other direction changes. Our calculation is two-dimensional and 
directly describes confocal arrays of cylindrical lenslets that are invariant to 
translation in one direction (for example the y direction; Fig. [2^a)), but the 
calculation also describes — separately — the relevant lateral projection^*"] 
of light passing through rectangular arrays of lenslets with rectangular aper- 
tures (Fig. (2^b)). If Q xz is the fraction of the power of the incident light that 
undergoes the correct direction change in the xz projection, and ( yz is the 
corresponding fraction in the yz projection, then the overall fraction of the 
power of incident light that undergoes the correct direction change is ( xz ( yz . 
Our calculation makes approximations by assuming that there is no dead area 
between neighbouring lenslets, and that each lenslet redirects light rays like 
an ideal thin len^" 2 "] Our results are important as they can be used to answer 
the question whether or not CLAs, along with realising previously forbidden 
refraction qualitatively, can also overcome previously established quantitative 
limits. 



If the lenslet arrays are parallel to the xy plane and periodic in the x and y 
directions, then the relevant lateral projections are into the xz and yz planes. 
2 This assumption implies that the individual lenslets have flat fields. 



3 



2 Confocal lenslet arrays (CLAs) 

CLAs (Fig. [T]) are formed by two parallel arrays of lenslets (or microlenses) , 
separated by the sum of their focal lengths [6j. In the simplest case, which 
we consider here, the lenslets' optical axes are perpendicular to both array 
planes, and each lenslet in one array has a corresponding lenslet in the other 
array with which it shares an optical axis. Corresponding lenslets can be seen 
as tiny two-lens telescopes, so CLAs are simply arrays of telescopes. Note that 
there is a generalisation of CLAs [7], but we do not consider it here. 

When a light ray is transmitted through a telescope, its transverse positions on 
entering and on exiting the telescope are generally different. This transverse 
offset occurs also in CLAs, but by miniaturising the telescopes the offset can 
be made smalQ Each telescope then acts like a pixel of the window formed by 
the CLAs; under the right conditions, the pixellation can be as unnoticeable 
as a computer monitor's. That this approach works has been demonstrated 
experimentally [13J. 

The change in direction of light rays that pass through a telescope consisting of 
lenses with focal lengths fi and fi can be described by the following equation, 
which describes the relationship between the angles of the light ray with the 
optical axis on the two sides of the telescope, a and (3 (Fig. [I]): 

fi tana = — / 2 tan (3. (1) 

As transmission through CLAs is transmission through telescopes, CLAs there- 
fore refract light rays according to this law of refraction, which is remarkably 
similar to Snell's law [BJ. 

The law of refraction given by Eqn ([!]) is interesting as it leads to perfect 
imaging (Fig. [ljb)) [6j. Obviously, if it is realised experimentally with CLAs, 
then the imaging can only be as good as the offset is small, and as good as 
the lenslets redirect light like ideal thin lenses. The case f\ = f% corresponds 
to refraction between media with refractive indices of equal magnitude and 
opposite sign [T4"][T5] . 

Strictly speaking, Eqn applies only to light rays that pass through cor- 
responding lenslets, that is, light rays that exit the same telescope that they 
entered [IE] . Provided this is the case, the direction change is independent of 
the precise position where the ray hits the first lenslet. The direction change 
such light rays undergo is called standard refraction; light rays that enter one 
lenslet and exit a non-corresponding lenslet (like the dashed light ray in Fig. 



The limit of useful miniaturisation has been exceeded when wave-optical effects 
begin to dominate. 
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[lja)) undergo a different direction change called non-standard refraction [T6] . 

In this paper we consider only CLAs with a particularly simple geometry: 
the aperture width of each lenslet is the same in both arrays, and there is no 
dead area between neighbouring lenslets (i.e. the centre-to-centre separation 
between neighbouring lenslets equald the aperture width). We describe such 
CLAs in terms of the dimensionless focal-length ratio [UJ, 

v = -T> W 
ji 



and the f- number of the left lenslet, 



N = A (3 ) 
where r is half the lenslets' aperture width. 

Our analysis considers light travelling from left to right, and assumes that the 
effect of transmission through the two lenslets can be described by first taking 
into account the effect of the lenslet with focal length f\ and then that of the 
lenslet with focal length f'2. In other words, it assumes that the lenslet with 
focal length f\ is to the left of the lenslet with focal length f'2, or that the two 
are in (or imaged into) the same plane, which is only the case if fi + $2 > 0. 
In terms of 77 and N, this becomes 

— 77) > (4) 

(as r > 0, by definition). The values of 77 and N are therefore restricted to 
combinations that satisfy this condition. 



3 Light undergoing standard refraction 



First we calculate the fraction of incident power that undergoes standard 
refraction. For simplicity, we consider here (and in the following sections) the 
incident light to consist of parallel rays, all of the same brightness; the angle 
of incidence is a. Wave-optically, this is a uniform plane wave. 



3.1 The case \q\ < 1 



First we treat the case \r}\ < 1, which corresponds to f\ > l^l- For this case, 
condition Q is satisfied provided that AT > 0. It is convenient to treat the 
ranges —1 < rj < and < 77 < 1 separately. 
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Fig. 3. (a) Geometry of and (b) critical beam paths through confocal lenslet arrays 
for the parameter range — 1 < r] < 0. (a) A collimated light beam enters a pair 
of corresponding lenslets from the left, inclined by an angle a with respect to the 
common optical axis. The longitudinal cross-sectional area of the beam is shown 
in solid light red; its edges are marked by solid red lines. The half-width of the 
lenslets is r. Counter-clockwise angles are positive; upwards and rightwards distances 
are positive, (b) Longitudinal cross-sectional areas for beams incident at critical 
angles. The part of the beam that passes through corresponding lenslets is filled in 
solid light red and its edges are marked by solid red lines. The part of the beam 
that enters through the left lenslet and exits through the lenslet above (below) the 
corresponding right lenslet is filled with green horizontal stripes (turquoise polka 
dots). 



Before calculating the fraction of the light that passes through the first (left) 
lenslet, which then passes through the corresponding (right) lenslet, we briefly 
discuss what we expect to happen for different values of the angle of incidence. 
Key throughout is the dependence of the light in the plane of the right lenslet 
array, described by the parameters a and b, on the angle of incidence, a (see 
Figs[3^a) and[4^a)). To derive this relationship, we start with the equations 



tan 7 



r + x 

~7T 



a — x 



tan (3 



r — x 



h 



b + x 



(5) 
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Fig. 4. (a) Geometry of and (b) critical beam paths through confocal lenslet arrays 
in the parameter range < r\ < 1. See the caption of Fig. [3] for further details. 

and so 



rf 2 + x(f 1 + f 2 ) 
h 

h 



(6) 
(7) 



With x = f\ tana, these expressions become the sought-for equations describ- 
ing the dependence of a and b on a: 



a = rf 2 / fi + (fi + f 2 ) tan a, 
b = rf 2 / f x - (/x + f 2 ) tan a. 



(9) 



For normal incidence (a = 0°), a and b take on the value a = b = rf 2 / fx- 
Specifically, \a\, \b\ < r, which means that all of the light that enters through 
the left lenslet exits again through the corresponding right lenslet. 

As a is increased to a "critical angle" a c2 , the upper edge of the beam hits the 
upper edge of the right lenslet. In Ref. [16], a c2 is called the "second critical 
angle of incidence"; we will encounter the other critical angles in due course. 
For — 1 < 7] < 0, a c2 is reached when a = r (see Fig. p[b)). Substitution into 
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Eqn (p| gives 



c2 _ K/i - h) _ i i + V 



For < ?7 < 1, a c2 is reached when b = — r (see Fig. |4](b)) ; substitution into 
Eqn @ gives 

tW 2 = -^ = ^ > 0). (11) 

Any further increase in a means that part of the beam misses the corre- 
sponding right lenslet. Another critical angle, a c3 (the third critical angle of 
incidence, according to the definition in Ref. [15]), is reached when the lower 
edge of the beam passes through the upper edge of the right lenslet. For 
— 1 < r] < 0, this happens when b = —r (see Fig. [3|b)); substitution into Eqn 
Q yields 

For < rj < 1, it happens when a = r (see Fig. [4](b) ) ; substitution into Eqn 
(|8| gives 

For any value a > a c3 , none of the beam passes through the corresponding 
right lenslet. 

It is perhaps worth noting that the expressions for a c3 for —1 < 77 < (Eqn 



(12)) and for a c2 for < rj < 1 (Eqn (11)) are identical, and so are the 



expressions for a for — 1 < rj < (Eqn (10)) and for a for < rj < 1 



(Eqn (13)). 



So far we have only discussed positive angles of incidence. What happens 
for negative values of a? For symmetry reasons, the same as what happens 
for positive values of a with the same modulus. Starting again from normal 
incidence and this time decreasing the value of a, the beam will again start to 
be clipped as another critical angle, a c ~ 2 , is reached. If a is decreased further, 
the fraction of the beam intersecting the corresponding right lenslet decreases 
until it reaches zero, at the critical angle a c ~ 3 . From the symmetry, and from 
Figs[3^b) and|i](b), it can be seen that 

a c ~ 2 = -a c2 , a c ~ 3 = -a c3 . (14) 



Now we calculate the power fraction of a uniform plane- wave beam of light that 
passes through corresponding lenslets. As before, the beam is incident from 
the left; we consider the part of the beam that has passed through a specific 
lenslet in the left lenslet array. In the plane of the corresponding right lenslet, 
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the height of the beam is |a + 6|. In the diagrams shown in Figs[3^a) and|4j(a), 
this is also the part of the corresponding right lenslet that is illuminated by 
this beam. The reason is that those diagrams are drawn for cases in which all 
light that passes through the left lenslet also passes through the right lenslet, 
that is, an angle of incidence a between a c ~ 2 and a c2 . For angles of incidence 
outside this range, the height of the second lenslet that is illuminated by the 
beam is less than the height of the beam, as some of the beam now misses 
the second lenslet. Generally, we can calculate the fraction £ of power that 
undergoes standard refraction as the fraction of the height of the right lenslet 
that is illuminated by the beam and the height of the beam in the plane of 
the right lenslet, \a + b\. 

For angles of incidence between 0° and the second critical angle, a c2 , the 
fraction ( is simply 1; above the third critical angle, a c3 , ( = 0. For angles of 
incidence between a c2 and a c3 , the top of the beam in the plane of the right 
lenslet is above the top edge of that right lenslet. The height of the beam 
illuminating the right lenslet is therefore r + b in the case — 1 < rj < and 
r — a in the case < rj < 1, and so 

C = L^ = V-l-2Nij,-l)^a { _^^ ol (15) 

s -a - b 2r] v - / - ; v j 

For negative angles of incidence, the fraction is the same as for a positive angle 
with the same modulus. These results can be summarised as follows: 



P 



if la I < a 



<2 



V+s g n( V )[l + 2N( V -l)t,n\a\] [{ ^2 < ^ < ^ ^ 

ifa c3 <|a|, 



where sgn(x) is the signum function, i.e. 



sgn(xj 



+1 if x > 0, 
if x = 0, 
-1 if x < 0. 



Fig. [5] shows an example plot of ( as a function of the angle of incidence, a. 
The role of the critical angles a c2 and a c3 as can clearly be seen: ( = 1 for 
— a c2 < a < a c2 , and ( = for |a| > a c3 . 

Fig. [6] shows a number of such curves, plotted for different values iV and 77. It 
is possible to observe compatibility with a number of trends in those graphs, 
for example the monotonic growth of a c2 between 77 = —1, where a c2 = 0, and 
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Fig. 5. Example of a plot of £ as a function of a. The graph is plotted for r\ = —0.5 
and N = 0.5 (e.g. f\ = 1, f% = 0.5, r = 1). The critical angles are indicated on the 
a axis; their values are a c2 = 18.4° and a c3 = 45°. 
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Fig. 6. Dependence of C( a ) on N, plotted for different values of r] in the range 
— 1 < r) < 1. The curves are plotted for N values between 0.25 and 2 in steps of 
0.25. 



f] = 0, where a c2 = arctan(l/2A r ) while a c3 stays constant for any fixed value 
of N; the equality of a c2 and a c3 at rj = 0; and the monotonic growth of a c3 
between 77 = and rj = 1 while a c2 remains constant. 
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Fig. 7. (a) Geometry of and (b) critical beam paths through confocal lenslet arrays 
in the parameter range r\ < — 1. See the caption of Fig. [3] for further details. 



Another — not particularly surprising, but nevertheless important — trend 
that can be seen in Fig. [6] is the narrowing of the ((a) curve as \N\ increases. 
This trend is easily explained: an increase in the modulus of the f-number 
corresponds to a decrease in the aperture width of all lenslets, and therefore 
a decrease in r, which makes it easier for light that has passed through one 
lenslet to miss the corresponding lenslet. This corresponds to a decrease in 
the field of view, which we take to be the range of angles between — a c3 and 
+a c3 , i.e. the angle range for which standard refraction occurs. 



3.2 The case \q\ > 1 



Next we consider the case \rj\ > 1, i.e. |/i| < j'2, which we treat analogously 
to the previous case (section j3T). Condition Q is now satisfied provided that 



N < if 77 > 1, or provided that N > if 77 < -1. 

The key quantities this time are a' and b', which describe the area of the left 
lenslet through which a light beam with angle of incidence a can enter for all 
of it to exit through the corresponding right lenslet (Figs[7^a) and[8^a)). We 



11 



(b) 

a c " 3 










7 


ill 





Fig. 8. (a) Geometry of and (b) critical beam paths through confocal lenslet arrays 
in the parameter range rj > 1. See the caption of Fig. [3] for further details. 

derive expressions for these quantities from the equations 



tan 7' 



b 1 + x 



r — x 



fi h 
and with x — f\ tan a (as before) we find 



tan/3' 



a — x 



fi 



r + x 



(19) 



a' = ^-[r + (/1 + / 2 )tana] 
h 

V = fy[r- (fi + h) tan a] 
h 



(20) 
(21) 



Starting from normal incidence (a = 0°) and increasing the angle of incidence, 
what happens qualitatively? For normal incidence, a' — b' — rf\j f^- In analogy 
to the previous case, |a'|,|6'| < r, but unlike in the previous case this now 
implies that not all the light that enters through the left lenslet exits through 
the corresponding right lenslet: only the light rays that enter the left lenslet 
in the shaded region (which, for a = 0°, is centred on the optical axis) exit 
through the corresponding right lenslet. As a is increased, the part of the 
left lenslet through which light rays have to enter to then exit through the 
corresponding lenslet moves upwards, but does not change in size until one 
of its edges (the upper side in Fig. H(b), the lower side in Fig. |8^b)) reaches 
the edge of the lenslet. The angle for which it happens is defined as the first 



12 



critical angle of incidence, a cl [16]. If rj < — 1, it can be calculated from the 
condition a' = r (see Fig. [7(b)), giving 



,1 r fi- f 2 1 1 + rj 



If T) > 1, the relevant condition is a' = — r (see Fig. [Sj^b)), and so 

t & na cl = -j = -~ (v >!)■ (23) 



As the angle of incidence is increased further, the part of the lenslet through 
which light rays have to enter to exit through the corresponding lenslet de- 
creases in size until its size is zero. If rj < —1, this happens when V = —r (see 
Fig. |7[b)), i.e. when 

tana c3 = -^ = ^ < -1); (24) 

if T) > 1, the condition is b' = r (see Fig. [81(b)) , and so 



c3 r fx - h 1 1 + rj 



^—.tvfr™^ (25) 



As before, negative values of a behave like their positive counterparts with 
the same modulus. The relevant critical angles (see Figs [7(b) andgb)) are 

oT x = -a c \ a c ~ 3 = -a c3 . (26) 



The power fraction of a uniform plane-wave beam that passes through corre- 
sponding lenslets is now the width of the left lenslet through which the beam 
passes, divided by the width of the lenslet. For a c_1 < a < a cl , this fraction 
is | a' + b'\/(2r) = I/1//2I = M -1 ; for a < a c ~ 3 and a > a c3 , it is zero. In 
between, i.e. for a c ~ 3 < a < a c ~ x and a cl < a < a c3 , it is given by 



C= ^h 2^ ^ H = 2, ^ * " X) 



t /2-/1 , + , , r ? + l + iV(r / -l)tan|a| 



(27) 
(28) 
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Fig. 9. Example of a plot of £ as a function of a for the case |r/| > 1, here rj = 2 and 
iV = —0.5. The critical angles are a cl = 45° and a c3 = 71.6°. 

In summary, 



\m 



-i 



if Id < a cl 



c = 



T?+sgn(T ? )[l+JV(7 ? -l) tan |q|] •£ c l < I I < a c3 
2t? — II — 

if a 63 < Id. 



(29) 



Fig. [9] shows an example of a plot of ( as a function of a, calculated according 
to Eqn ( p9[ ). The main difference with the corresponding plot for the case 
1^1 < 1 (Fig. [5j is that the maximum value C reaches is not 1, but l/\f]\- As 
before, the role of the critical angles, this time a cl and a c3 is apparent in the 
graph. 



Fig. [10] shows a number of such curves for different values N and rj, grouped 
together by their value of rj. Obviously, the maximum value of these curves 
is the same for values of i] with opposite signs, and it approaches the value 1 
as r] approaches ±1. The decrease in the field of view, i.e. the range of angles 
between — a c3 and +a c3 , with increasing \N\ is again apparent. 



4 Light undergoing non-standard refraction 



We now turn our attention to non-standard refraction |16j . Non-standard re- 
fraction is defined as the direction change of light that passes through non- 
corresponding lenslets, i.e. lenslets that do not share an optical axis. 

First we derive the law of refraction for non-standard refraction. We can do this 
by treating it a special case of refraction due to generalised CLAs |7] : arrays of 
lenslets in which the telescopes formed by pairs of corresponding lenslets have 
been modified, but always such that they continue to share a common focal 
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Fig. 10. Dependence of C(°0 on N, plotted for different values of rj in the range 
\rj\ > 1. The curves are plotted for N values between 0.25 and 2 for q < — 1, and 
between —0.25 and —2 for rj > 1, in steps of 0.25 in both cases. 



plane. Suitable modifications can include, for example, a sideways translation 
of one of the lenslets. In Ref. [7], such a translation is defined in terms of 
a dimensionless parameter 5 = d/ ft, where d is the translation distance. If 
the sideways translation is the only generalisation of standard CLAs (Fig. 
[TJ, then the law of refraction describing transmission through the CLAs is 
tana = 5 + rj tan (5. As expected, for S = 0, this equation becomes the law of 
refraction for standard refraction, Eqn (11]). As a non-corresponding lenslet can 
simply be seen as the corresponding lenslet, but translated sideways through a 
distance d = m2r, where m is an integer, non-standard refraction can be seen 
as a special case of refraction with generalised CLAs with 8 = m2r/fi = m/N 
(and 6 = 0, i.e. no rotation). Therefore the law of refraction for non-standard 
refraction through non-corresponding lenslets is 



tana 



m 
N 



+ T/tan/3. 



(30) 



Figs 11 and 12 show graphs of the angle of refraction, /3, as a function of the 
angle of incidence, a, for a few combinations of rj and m/N. 
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Fig. 11. Law of non-standard refraction, Eqn (30), plotted for various values of 77. 
The curves are shown only in the range of incidence angles for which a non-zero 
fraction of the light exits through non-corresponding lenslets shifted by m lenslets, 
i.e. that between the critical angles a^" 3 and c£?. In all cases, N = 1; for plots of 



the dependence on N, see Fig. 12 



Next, we derive the power fraction of incident plane-wave light that undergoes 
non-standard refraction with a given value of m. We refer to Fig. 13, which 
sketches the geometry of light undergoing non-standard refraction with n = 1 
in the four cases (—1 <?7<0, < 77 < 1, 77 < —1, and 77 > 1) to which Figs[3l 



|4j [7] and [8] refer. Like in the calculation of the power fraction that undergoes 
standard refraction (section [3]) , the quantities a, b, a' and b' (see Fig. 13) are 
key to our current calculation, a and b are still given by Eqns (JSh , d9|), but a' 
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and V are now dependent on m, which we indicate with a subscript m: 




For m — 0, these expressions reduce to the ones we used in section 3.2[ Eqns 
(20) and (21), as one would expect. The calculations in this section generally 
follow the same outlines as those in sections |3.1| and |3.2[ so we keep them 
brief. 



The diagrams in Fig. 13 are drawn for values of the angle of incidence, a, 
that have been chosen such that the maximum fraction of the incident light 
undergoes non-standard refraction with m — 1. For \rj\ < 1, this power fraction 
is 1; for \rj\ > 1, it is as before. 



As a is increased, the transmitted power fraction stays constant until the beam 
hits the top edge of the right lenslet (for \rj\ < 1) or the left lenslet (for \rj\ > 1). 
This happens when a = (2m + l)r if — 1 < 77 < 0; when b = —(2m + l)r if 
< T] < 1; when a' m = r if 77 < — 1; and when a' m = — r if r\ > 1. We take 
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Fig. 13. Geometry of light undergoing non-standard refraction. In the examples 
sketched here, light enters the left lenslet but does not exit through the correspond- 
ing lenslet on the same optical axis. Instead, it exits through the mth lenslet above 
the corresponding lenslet. The diagrams are drawn for m = 1. 

these conditions to define critical angles; for \r]\ < 1, 



tana< c2 ' m 



1^2^+1+2 if _ i < ^ < o, 



2N 1-r) 
1 2m+l~r) 



, 2N 1— r? 11 U — — X > 



for \r)\ > 1, 



(33) 



tana cl ' m = I 2N ^ ' 34) 
2N i-n I — 



As a is increased further, the power fraction reduces until it reaches zero, 
when b = — (2m + l)r if — 1 < rj < 0; a = (2m + l)r if < 77 < 1; b' m — — r if 



18 



f] < — 1; b' m = r if 77 > 1. These again define critical angles: 



tana 



c3,m 



1 


2m+l-»7 


27V 




1 


2m+l+T) 


27V 


1-7, 


1 


2m+l— »j 


27V 




1 


2m+l+rj 


I 27V 


1-1 



if - 1 < 77 < 
if < 77 < 1, 
if 77 < -1, 
if 77 > 1. 



(35) 



Starting again from incidence angles for which the maximum fraction of the 
incident light undergoes non-standard refraction (angles such as the ones for 



which Fig. 13 is drawn), we now briefly consider what happens when the angle 
of incidence is decreased. 

The power fraction stays constant until the beam hits the bottom edge of the 
right lenslet (for |?7| < 1; see Figs [31(b) and|4^b)) or the top edge of the left 
lenslet (for |?7| > 1; see Figs[7^b) and[8^b)). This happens when b = — (2m — l)r 
if — 1 < 77 < 0; when a = (2m — l)r if < 77 < 1; when b' m = r if rj < —1; and 
when b' m = — r if 77 > 1. For \r]\ < 1, we call the angle for which this happens 



-' i,m , and calculate it to be 



a 



c—2,m 



1 2m-l-i) -r 
27V 1-r) 11 



1 < 77 < 0, 



i^sm-i+a if < 77 < 1. 



(36) 



27V 1-n 



As 2m - 1 = 2(m - 1) + 1, a 
the relevant angle a 



c—2,m 



a: 



c3,m— 1 



in both cases. For I77I > 1, we call 



c— l.m 



and we find that 



a 



-l.m 



^?^ifr7<-l, 



2N 

1 2m+l-t) 
27V 1-n 



if 77 > 1. 



(37) 



As the angle of incidence is decreased further, the power fraction decreases 
until it reaches zero, when a = (2m — l)r for — 1 < 77 < 0; b = —(2m — l)r for 
< 77 < 1; a' m — — r for 77 < — 1; a' m — r for 77 > 1. The corresponding critical 
angle, a c ~ 3 ' m , is then 



a 



c—3,m 



1 


2m- 


14 


-v 


27V 


1- 


-v 




1 


2m- 


1- 


V 


27V 


1- 


-n 




1 


2m- 


14 


-v 


27V 


1- 


->i 




1 


2m- 


1- 


V 


27V 


1- 


->i 





if - 1 < 77 < 0, 
if < 77 < 1, 



(3* 



Finally, we are ready to calculate the power fraction of a uniform plane wave 
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incident at an angle a that enters through a lenslet in the left array and exits 
through the lenslet in the right array that is m lenslets above that correspond- 
ing to the entrance lenslet. We call this power fraction ( m . 

We start with the easy cases. For angles of incidence, a, below a^ 3 and 
above a^, Cm = irrespective of the value of 77. If \r)\ < 1, then Cm = 1 for 

2 < a < a m! if \V\ — 1; t nen Cm = l/\v\ f° r a m 1 < a < 

The remaining cases are slightly more complicated, and there are a number of 
them. They can be calculated from the following expressions. 

(1) For -1 < T) < 0, Cm = [(2m + l)r + b]/(a + b) if a£? < a < and 
C m = [-(2m - l)r + a] /(a + 6) if a^ 3 < a < a c ~ 2 ; 

(2) for < T) < 1, Cm = [(2m + l)r - a]/(-a - 6) if a% < a < a c3 and 
C m = [-(2m - l)r - 6]/(-a - 6) if a4" 3 < a < <~ 2 ; 

(3) for V < -1, Cm = [r + VJ/(2r) if < a < < and Cm = [r + a' m ]/(2r) 
if a^; 3 < a < a c ~ 2 \ 

(4) and for V > 1, ( m = [r - b' m \/{2r) if a% < a < a% and Cm = [r - <J/(2r) 
if a c ~ 3 <a< a c m 2 . 



The results can be summarised as follows. For — 1 < r] < 1, 



r?+sgn(r?)[-l+2m-2(l-q)Af tan a] 
1 

r?+sgn(r?)[-l-2m+2(l-»?)Af tan a] 




if a < a 



c-3 



if a 
if a 



c-3 



c-2 



< a < a 

< a < a 



c-2 
m ) 

c2 
m ' 



if o-t < a < ot%, 
if < «; 



(39) 



for \r)\ > 1, 





?7+sgn(q)[l-2m+2(l-r?)jVtan 



if a < a 



c-3 



if oi c 3 < a- < o' 1 

2r\ m — " — "m 5 

J_ 



if a m 1 < a < a m , 



t ? +sgn(? ? )[l+2m-2(l-T ? )jVtana] -r c l < < c3 
2ri m — " — "m; 

if aS 3 < a. 



(40) 



Fig. 14 shows plots of Cm as a function of a for various different values of 77. A 
number of things can be seen in those graphs. Perhaps most notably the curves 
suggest that the values of C for all values of m add up to one for all values 
of 77 and all angles of incidence a between —90° and +90°, and indeed we 
confirmed this numerically for a number of cases. It can also be seen that the 
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Fig. 14. Power fraction £ m (Eqns (39) and (40)), plotted as a function of the angle 
of incidence, a, for m = 0, ±1, ±2 (respectively represented by lines that are solid, 
dashed, and dotted) for various values of r\. All curves were calculated for \N\ = 1, 
i.e. N = —1 for rj = 2 and N = +1 for all other cases. 

peak in the ( m curve moves to greater values of a as m is increased (provided 
the other parameters are unchanged). 



All the curves in Fig. 14 were plotted for \N\ 



1. Fig. 15 shows how the 
curves change as only is varied. As \N\ is increased, the peaks become 
narrower, and the peaks corresponding to greater values of |m| move inwards. 



5 Similarities with multiple refraction and reflection of cold atoms 
in light fields 



It is interesting to note that the law of refraction for light passing through 
non-corresponding lenslets, Eqn (30), allows a link to recent work on cold 
atoms. This works as follows. 

With a suitable choice of light fields, the center-of-mass motion of colds atoms 
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Fig. 15. Plots of the power fraction £ m (Eqns (39) and (|40|) as a function of the 
angle of incidence, a, as N is varied. Graphs for m = (solid lines), m = ±1 
(dashed lines) and m = ±2 (dotted lines) are shown. All graphs were calculated for 
ri = -0.5. 

can be described by the non-Abelian vector potential term proportional to 
the spin- 1/2 operator [TTfT 8 ,19,20] . The atomic motion is then characterised 
by two dispersion branches containing the areas of both positive and nega- 
tive slopes. In the low-energy part of the spectrum for each frequency, there 
are two wave vectors corresponding to the positive and negative group veloc- 
ity, respectively. As a consequence, the incident atomic wave may split into 
two reflected waves at a barrier, one that undergoes specular reflection, and 
an additional wave that undergoes non-specular reflection [20]. Note that a 
similar kind of double reflection can occur also for electrons affected by the 
Rashba spin-orbit coupling [2"T1I2"2"] . On the other hand, the transmitted wave 
experiences the negative refraction for small angles of incidence [19], whereas 
for large angles of incidence an additional positive refraction takes place [23] . 
The former negative refraction is similar to that taking place by passing the 
light through the corresponding lenses, whereas the latter ordinary refraction 
corresponds to the light passing through non-corresponding lenslets. 



6 Conclusions 



We have calculated the fraction of optical power that undergoes standard 
refraction on transmission through CLAs. The calculation requires different 
cases to be treated separately, making it slightly cumbersome. The fruit of 
this drawn-out labour is a number of equations describing the transmission of 
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optical power through CLAs, which can be used for assessing the suitability 
of CLAs for potential applications. 

But the fraction of power transmitted through CLAs is only one of a num- 
ber of considerations in potential applications. An example of an additional 
consideration is the accuracy with which CLAs redirect light rays, which is 
determined by how close the effect of each individual lenslet is to that of an 
ideal thin lens. As lenslet arrays are usually arrays of very basic lenses (typi- 
cally bumps of glass or plastic), they work best for rays passing close to the 
centre, i.e. for high f-number. Our results show (not very surprisingly) that a 
CLA's field of view is lower for higher f-numbers, and so a compromise needs 
to be found that suits a particular application. Despite these complications, 
we believe that our results will be key to investigating the suitability of CLAs 
for different potential applications, specifically solar concentrators. 
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